Search for: All records

This paper presents a variational multiscale (VMS) based finite element method where the stabilization parameter is computed dynamically. The current dynamic procedure takes in a general structure/form of the stabilization parameter with unknown coefficients and computes them dynamically in a local fashion resulting in a dynamic VMS‐based finite element method. Thus, a static stabilization parameter with pre‐defined coefficients is not needed. A variational Germano identity (VGI) based local procedure suitable for unstructured meshes is developed to perform the dynamic computation in a local fashion. The local VGI based procedure is applied for each interior vertex in the mesh and unknown coefficients are first determined locally at each vertex, and subsequently, for each element a maximum value is taken over the vertices of the element. To make the current procedure practical, a coarser secondary solution is constructed from the primary coarse‐scale solution, which is done locally over a patch of elements around each interior vertex. Further, averaging steps are employed to make the local dynamic procedure robust. Currently, the new dynamic VMS formulation is applied to steady problems governed by the advection‐diffusion and incompressible Navier‐Stokes equations in both 1D and 2D to demonstrate its efficacy and effectiveness.

Numerical simulations of cardiovascular mass transport pose significant challenges due to the wide range of Péclet numbers and backflow at Neumann boundaries. In this paper we present and discuss several numerical tools to address these challenges in the context of a stabilized finite element computational framework. To overcome numerical instabilities when backflow occurs at Neumann boundaries, we propose an approach based on the prescription of the total flux. In addition, we introduce a “consistent flux” outflow boundary condition and demonstrate its superior performance over the traditional zero diffusive flux boundary condition. Lastly, we discuss discontinuity capturing (DC) stabilization techniques to address the well‐known oscillatory behavior of the solution near the concentration front in advection‐dominated flows. We present numerical examples in both idealized and patient‐specific geometries to demonstrate the efficacy of the proposed procedures. The three contributions discussed in this paper successfully address commonly found challenges when simulating mass transport processes in cardiovascular flows.